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We report transport measurements of gate-tunable Andreev bound states in a carbon nanotube 
quantum dot coupled to two superconducting leads. In particular, we observe clear features of two 
types of Kondo ridges, which can be understood in terms of the interplay between the Kondo effect 
and superconductivity In the first type (type I), the coupling is strong and the Kondo effect is 
dominant. Levels of the Andreev bound states display anti-crossing in the middle of the ridge. On 
the other hand, crossing of the two Andreev bound states is shown in the second type (type II) 
together with the O-tt transition of the Josephson junction. Our scenario is well understood in terms 
of only a single dimensionless parameter, /csT^*^/A, where T^*^ and A are the minimum Kondo 
temperature of a ridge and the superconducting order parameter, respectively. Our observation is 
consistent with measurements of the critical current, and is supported by numerical renormalization 
group calculations. 

PACS numbers: 73.23.-b, 73.63.Kv, 03.65.Yz 



In an Andreev reflection process at the interface be- 
tween a normal metal (or any other non- superconduct- 
ing region) and a superconductor (S), an electron in the 
normal region is converted to a Cooper pair necessitat- 
ing a reflection of a hole. Multiple Andreev reflections 
(MAR) play a central role in finite-bias transport through 
a non-superconducting region sandwiched by two super- 
conducting leads [IHS]. Contrary to MAR peaks at fi- 
nite bias, Andreev bound states (ABS) are formed as a 
result of coherent superposition of all possible Andreev 
reflection processes (to the infinite number). So far ABS 
have been observed either in equilibrium across two su- 
perconducting electrodes [6] , or in a system with only one 
superconducting lead [7l[8]. The interplay of ABS with 
the Kondo effect in S-quantum dot (QD)-S is an inter- 
esting issue, particularly in the context of the O-tt transi- 
tion [9HT4] . The "tt state" , exhibiting the reversed sign of 
the Josephson current of S-QD-S, originates from an un- 
paired electron spin in a strongly interacting QD [9l [10] . 
The usual "0 state" is recovered due to the Kondo screen- 
ing in the strong coupling limit. Theoretically, it has 
been shown that this O-tt transition is followed by the 
level crossing of two ABS [15]. However, simultaneous 
observation of ABS level crossing and the O-tt transition 
has never been achieved experimentally. 

In this Letter, we report a clear signature of a gate- 
tunable ABS in non-equilibrium transport through a car- 
bon nanotube QD coupled to two superconducting leads. 
In particular, we show that the interplay between ABS 
and the Kondo correlation plays a major role in trans- 
port. This leads to the two different prototypes of the 



Kondo ridges depending on the ratio ksT^''' / /\. T^^^ 
and A represent the minimum Kondo temperature of 
the ridge and the superconducting order parameter, re- 
spectively. "type-I" Kondo ridge with a stronger cou- 
pling {ksT^'^^ / A ^ 0.8) displays an anti-crossing of 
ABS, where the Josephson junction is always in the "0 
state". On the other hand, "type-II" Kondo ridge with 
a weaker coupling {ksT^'^^ /A ^ 0.8) shows a crossing 
of ABS, which is directly related to the O-tt transition of 
the Kondo-correlated Josephson junction. Our transport 
data clearly show two types of the Kondo ridge in a single 
sample. This scenario is confirmed by our measurement 
of gate-dependent Josephson critical current Ic along the 
ridge. In addition, our experimental result is supported 
by numerical renormalization group (NRG) calculations. 

Experimental setup. - Carbon nanotubes (CNTs) are 
grown by the conventional chemical vapor deposition 
method on Si02/Si (500 nm/ 500 /im) wafer [16 . A 
CNT is located by atomic force microscopy. Contact elec- 
trodes (10 nm/80 nm Ti/Al bilayer) are subsequently re- 
alized by electron-beam lithography and successive elec- 
tron beam evaporation processes [17J. Ti layer is used as 
an adhesion layer between the superconducting Al layer 
and CNT. The superconducting transition temperature 
Tc of the bilayer is about 1.1 K. All measurements are 
performed at a base temperature below 100 mK in a di- 
lution refrigerator. Two terminal DC measurements are 
done both in the current and in the voltage bias modes 
by using DC volt age /current source (Yokogawa GS200), 
as shown in Fig. 1(a) [18 . Al superconducting elec- 
trodes are switched to normal state by applying an ex- 
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FIG. 1: (a) Schematic of the measurement configuration. 
Highly doped Si wafer and Si02 layer are used as a back gate 
and an insulating barrier, respectively. Channel length of the 
CNT is designed to be 300 nm. Diameter of the CNT -1.5 
nm is measured by AFM. Differential conductance dl/dVsd 
plot as a function of Vsd and Vg is displayed (b) for normal 
(with magnetic field B ~ 1 kG), and (c) for superconducting 
states at zero field, respectively. The letters 'e' and 'o' denote 
the even and odd number states, respectively. Bias voltages 
corresponding to A/e and 2A/e (with A — 140 /xeV) are in- 
dicated by arrows in (c). Red curves in (b) correspond to 
dl/dVsd vs. Vsd in the middle of the Kondo ridges A - E. 



ternal magnetic field of ~ 1 kG at the base temperature. 

The normal state differential conductance dl/dVsd is 
plotted in Fig. 1(b), providing a typical even-odd behav- 
ior with pronounced Kondo peaks at the zero bias of odd- 
number ridges. From Fig. 1(b), we extract a charging en- 
ergy /7 ~ 3 meV and a level separation 5E ~ 3 meV. The 
Kondo temperature T^, estimated from the half width 
half maximum of the Kondo peaks [19], depends on the 
gate voltage. The minimum Kondo temperature T^*^ in 
the middle of each Kondo ridge is estimated to be 0.94, 
1.6, 2.6, 2.4 and 0.86 K for the Kondo ridges A to E, 
respectively. The tunnel coupling r(= F^^ + V^) (with 
tunnel coupling to the right (left) electrode Fi^(ri;,)) is 
estimated by fitting Tk data to the formula ksT^ = 
{UT/2)^/^exp[-7r[\4s^ -U^\]/8UT], where s is the energy 
level in QD tuned by gate voltages [20l |2T]. The esti- 
mated values of F range from 0.35 ~ 0.96 meV for ridges 
A-E. The asymmetry ratio j{= Tr/Tl) is obtained from 
the relation Gmax = {2e^/h) • 4VrVl/{Tr + F^)' ^22^, 
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FIG. 2: (a) Differential conductance plot for the two types of 
the Kondo ridge (Magnified view of Fig. 1(c) for ridge D and 
E). Red color represents the region of the negative differential 
conductance. Blue dashed lines {Vabs) are associated with 
the ABS-assisted transport (see the text), (b) Schematic di- 
agram of the major transport process forming the Andreev 
bound states at finite bias. This mechanism is expected to be 
pronounced in highly asymmetric barriers, where the Andreev 
bound states are pinned to the lead with stronger coupling, 
(c) The ABS level position Sa obtained from Eq. (1) is plotted 
as a function of Vg. (d) UbTk / ^ vs. Vg for ridges D (type-I) 
and E (type-II) , respectively. Filled circles and solid lines cor- 
respond to the experimental data and the theoretical fit with 
the formula ksTK = (UT /2Y^'^e^^[-Ti[\4e^ - U'^\]/WT], and 
[//A = 21.4. The dashed line at ksTK/A = 0.788 is the 
theoretically expected 'O-tt transition' line [13]. (e) Ic vs. Vg 
for ridges D and E, respectively. Solid lines are for an eye- 
guide. Ic-Vg curves for ridge E show sharp drop around the 
O-TT transition point denoted by the two vertical dot-dashed 
lines. 



where Gmax is the zero-bias linear conductance obtained 
in the middle of each Kondo ridge (not shown here) ^3] . 
The estimated 7 is 2.7, 16, 2.3, 12, and 3.1 for ridges A - 
E, respectively. This asymmetry plays an important role 
in ABS-assisted transport, as we discuss below. 

Main features of the finite-bias transport. - Differen- 
tial conductance with the superconducting leads is dis- 
played in Fig. 1(c). While the even number state shows 
the gate-independent (elastic) quasiparticle cotunneling 
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FIG. 3: NRG results for the two types of the Kondo ridge: 
(a) ABS level position £a (in unit of A) and (b) the Joseph- 
son critical current Ic as a function of the QD energy level 
£d- Both in (a) and (b), coupling constants F/A = 4.9 (left 
panels) and F/A = 3.0 (right panels) are used, which corre- 
spond to the experimental values of the ridges D (type-I) and 
E (type-II), respectively. Experimentally extracted value of 
[//A= 21.4 is used in both cases. /csT^'^'/A are 1.5 and 0.5 
for ridge D and E, respectively. 



{^ysd — ±2A) [24] together with a weak single Andreev 
reflection (eVsd = peaks, the odd-number state dis- 
plays a rich subgap structure (|eVs(i| < 2A), which origi- 
nates from the interplay between the Kondo effect and su- 
perconductivity [19 . The most prominent feature is the 
two different prototypes of the Kondo ridges. This can be 
seen more clearly in Fig. 2(a), a magnified view of ridges 
D and E in Fig. 1(c). Fig. 2(a) displays a strong sub- 
gap transport. We focus on the main peaks {Vsd = Vabs) 
which vary from IcVabs'I = A to IgVabs'I = 2 A as a func- 
tion of the gate voltage (blue dashed lines). Such strong 
peaks cannot be understood in terms of the perturbative 
MAR peaks. Notably, this gate-dependent peak in the 
Kondo ridge evolves into the elastic quasi-particle cotun- 
neling peak at \eVsd\ = 2 A in the even valley. 

We attribute these main peaks to ABS-assisted trans- 
port, which is illustrated in Fig. 2(b). In highly asym- 
metric junctions, ABS are formed mainly between QD 
and a lead with lower barrier (left lead in Fig. 2(b)). The 
other lead with higher barrier would play a role of prob- 
ing ABS. In this picture, the gate-dependent bright peaks 
at Vsd = Vabs result from the alignment of ABS and the 
gap edge of the "probe" lead. Together with the fact 
that ABS is formed always in pair with the electron-hole 
symmetry [15], this gives the following relation between 
ABS energy (sa) and the peak position (Vabs)' 
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The gate-dependent ABS obtained from this relation 
is plotted in Fig. 2(c). ABS displays two distinct fea- 



tures depending on the type of the Kondo ridges. Type- 
I ridge (ridge D) with a stronger coupling displays an 
anti-crossing in the middle of the ridge where the Kondo 
temperature has its minimum value. On the other hand, 
type-II ridge (ridge E) shows two clear crossing points 
of ABS, which is related to the O-tt transition. The two 
different types are determined by the ratio /c^T^^^/A. 
In type-I ridge, this ratio is always larger than the crit- 
ical value {kBT^'^'/A 0.8) (left panel of Fig. 2(d)). 
In contrast, the Kondo temperature in type-II ridge has 
two transition points corresponding to O-tt transition in 
the right panel of Fig. 2(d)). 

Critical currents - Our scenario is further supported 
by the behavior of the Josephson critical current, Ic, 
(Fig. 2(e)) [25]. High Josephson current flows due to 
the Kondo-assisted transmission in the strong coupling 
ridge (type I), and the system is always in "0 state" (left 
panel of Fig. 2(e)). In type-II ridge, O-tt transition leads 
to a dramatic change in the behavior of Ic (right panel 
of Fig. 2(e)). As the Kondo effect is suppressed in the 
middle of the ridge, "tt state" appears, and it leads to 
a strong suppression of Ic- All these features are con- 
sistent with the behavior of the gate-dependent ABS of 
Fig. 2(c). This behavior of the supercurrent has also 
been reported in Ref. [21], without delving it into ABS 
in finite-bias transport. 

Numerical results - Our main observation in Fig. 2 is 
nicely supported by NRG calculation. Fig. 3(a) shows 
the calculated ABS level e^/A as a function of the QD en- 
ergy level Ed' It clearly displays two different prototypes, 
i.e. anti-crossing (left) and crossing (right) of ABS for pa- 
rameter values of F/A = 4.9 and 3.0, respectively (these 
values correspond to the experimentally extracted values 
for ridge D and E, respectively). Theoretical values of Ic 
are determined by the amplitude of the current-phase re- 
lation at zero bias voltage. Fig. 3(b) shows the QD-level 
dependence of Ic for the two types of the ridge. As in the 
experimental result of Fig. 2(e), large Josephson current 
is assisted by the Kondo effect for type-I ridge. Suppres- 
sion of the supercurrent is clearly shown for type-II ridge 
at "tt state" (right panel of Fig. 3(b)). This sudden drop 
of the supercurrent takes place exactly at the O-tt tran- 
sition point. All the features in the calculations match 
perfectly with the experimental observation in Fig. 2 [26] . 

Universality of ABS-assisted transport - We note that 
the gate-independent peaks for even number state at eVsd 
= ±2A can also be understood within our framework of 
ABS-assisted transport. Of course, these peaks can be 
identified as the elastic cotunneling of quasiparticles [24] . 
On the other hand, ABS are pinned to £a = ±A in the 
off-resonance limit of the even valley. In this case, ABS- 
assisted transport also provides peaks at eVsd = ±2 A, 
according to Eq. (1). That is, ABS peaks evolve into 
the quasi-particle cotunneling peaks in the off-resonance 
limit. Therefore, the main features of transport can be 
understood in a very universal way, ranging from the 
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FIG. 4: / — Vsd curves measured in a current bias (red line), 
and voltage bias (blue line) modes, respectively, in the middle 
of the Kondo ridge D (Vg=-8.255 V). Two switching currents, 
Ic and /a, appear at zero voltage and at V ^ 0.2 mV, respec- 
tively. The slope of the supercurrent branch is mainly due to 
the resistance of low pass filters (^ 4 kQ). Upper-left inset: 
differential conductance G = dl/dVsd vs. Vsd in voltage bias 
mode. Lower-right inset: magnified view of the I-V curve 
around zero bias. 

Coulomb blockade (off-resonance), the intermediate cou- 
pling (tt state), to the strong Kondo limit (0 state). 

Negative dynamic conductance and fine structures - 
Another evidence of ABS-assisted transport is shown in 
the negative dynamic conductance (NDC) region in the 
voltage bias mode (blue line in Fig. 4). It can be more 
clearly seen in dl/dVsd (upper- left inset of Fig. 4 and 
the red-colored region in Fig. 2(a)). Hysteretic switching 
current above Ic is found in a current bias at / = /q (for- 
ward) or / = lar (backward) depending on the direction 
of the current sweep. NDC and hysteretic switching cur- 
rent can be understood as a general feature when Cooper 
pair tunnels through a resonant state in a QD [27H29]. 
For highly asymmetric barriers (as illustrated in Fig. 
2(b)), resonant states formed with the lower barrier lead 
can be probed by the "probe" lead. However, as tunnel- 
ing barriers become more symmetric, ABS levels are not 
well defined at finite bias and NDC is expected to disap- 
pear. Actually, NDC and hysteretic switching current at 
I = la and / = lar are observed only for ridge B and D 
whose asymmetry ratio is very large (7 = 16 and 7 = 12 
for ridges B and D, respectively). When one of the two 
tunneling barriers is significantly higher than the other, 
our measurement configuration is equivalent to scanning 
tunneling microscope (STM) measurement with a super- 
conducting tip. Thus, the dynamic conductance density 
plots for ridge B and D reflect the density of states in the 
CNT-QD as shown in Refs. [SHH]. 

Finally, we briefly discuss the fine structures in the 
differential conductance at lower voltage, \eVsd\ < A in 
Fig. 1(c) and Fig. 2(a). Although a quantitative anal- 
ysis for this is beyond our scope, we notice that the 



smaller peaks at \eVsd\ < A greatly resemble the gate 
dependence of the main peaks at Vsd = Vabs for both 
types of the Kondo ridge. Therefore, we speculate that 
these small peaks originate from the combination of ABS 
(mainly formed with the lower barrier contact) and single 
or multiple Andreev reflections (with the higher barrier 
contact). 

Conclusion - In conclusion, gate tunable ABS are re- 
ported in I-V measurement configuration in an Al-CNT- 
Al Josephson junction. The observed differential conduc- 
tance shows the two distinct types of the Kondo ridges as- 
sociated with ABS. ABS displays crossing (anti-crossing) 
behavior, which is the main characteristics of the O-tt 
transition (0 junction) tuned by a gate voltage applied 
to the QD. This feature is also consistent with a mea- 
surement of the gate-dependent critical current, and is 
confirmed by a numerical renormalization group calcula- 
tion. 
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